On the relation of quark confinement and chiral symmetry breaking 
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We study the phase diagram of two flavour QCD at imaginary chemical potential in the chiral limit. 
To this end we compute order parameters for chiral symmetry breaking and quark confinement. 
The interrelation of quark confinement and chiral symmetry breaking is analysed with a new order 
parameter for the confinement phase transition. We show that it is directly related to both, the 
quark density as well as the Polyakov loop expectation value. Our analytical and numerical results 
suggest a close relation between the chiral and the confinement phase transition. 
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Introduction - Quantum Chromodynamics (QCD) at 
finite temperature and density is a very active area of re- 
search. The equation of state of QCD and, in particular, 
the nature of the transition from the hadronic phase with 
broken chiral symmetry to the chirally symmetric decon- 
fined quark-gluon plasma phase is of great importance for 
a better understanding of the experimental data, e.g. 

For full QCD with dynamical quarks one expects the 
confinement phase transition to be a crossover as quarks 
explicitly break the underlying center symmetry of the 
gauge group. The nature of the chiral phase transition 
primarily depends on the value of the current quark mass, 
which explicitly breaks chiral symmetry, as well as on 
the strength of the chiral anomaly 0. While the con- 
finement phase transition is driven by gluodynamics, the 
chiral phase transition is governed by strong quark in- 
teractions. Hence, it is a highly non-trivial observation 
that both lie remarkably close at least for small quark 
chemical potentials 3 ■ An understanding of this in- 
terrelation is subject of an ongoing debate. 

The chiral properties of the theory are accessi- 
ble through the spectrum of the Dirac operator Q. 
Recently it has been shown that also the confining 
properties of the theory can be accessed through its 
spectrum 0, 0, H, In this Letter we analyse the 
deconfinement and chiral phase transition and their 
interrelation both analytically and numerically for two 
flavour QCD in the chiral limit. The interrelation is 
studied with the help of a new order parameter for the 
confinement phase transition which is related to the 
quark density. The numerical analysis is performed with 
functional renormalization group (RG) methods, for 
reviews see [U , and suggests that the deconfinement 
phase transition is indeed correlated with chiral symme- 
try restoration. 

Order parameters - The Polyakov loop variable L(x), 
L{S) = — trP(f) with V{x) = p e'9 I<f dt Aoit,x) ^ 

in QCD with Nc colors and infinitely heavy quarks is re- 
lated to the operator that generates a static quark p^ . 
In II]) the trace is evaluated in the fundamental repre- 



sentation, and P stands for path ordering. We can inter- 
pret the logarithm of the expectation value (L) as half 
of the free energy Fqq of a static quark-anti-quark pair 
at infinite distance. Moreover the expectation value (L) 
is an order par ameter for center symmetry of the gauge 
group, e.g. [Ij]- To see this, we consider gauge trans- 
formations Uz{t,x) with U~^{0, x)Uz{P, x) = z, where 
z G Z is an element of the center Z of the gauge group. 
Under such a transformation the Polyakov loop is multi- 
plied with a center element z, L(x) —^ zLix). Hence a 
center-symmetric confining disordered ground state with 
-Fq — > oo is ensured by (L) = 0. In turn, deconfinement 
with Fg < oo is signaled by {L) ^ 0. This implies center- 
symmetry breaking in the ordered phase. 

It follows immediately that any observable, which 
transforms non-trivially under center transformations, 
serves as an order parameter. This has been exploited 
in d, 0, [1, [§1 , where the spectral properties of the Dirac 
operator have been related to the expectation value of 
the Polyakov loop. The relation stems from the obser- 
vation that, in contradistinction to the gauge fields, the 
periodicity properties of the quark fields change under 
application of a gauge transformation C/z, 

^^'[t + (i,x) = -ziij^'[t,x) , with = (2) 

A straightforward generalisation of the boundary condi- 
tions of the quarks yields 



(3) 



This includes ([2]) with the center phases z = Ue^'^'^^, 
e.g. 9z = 0,1/2 in SU{2), and 9^ = 0,1/3,2/3 in SU{i). 
Quarks with the boundary conditions ([3|) can be rewrit- 
ten in terms of quarks with physical anti-periodic bound- 
ary conditions, 

,l;e{x)=e^''^'^l^i!{x) with i}{x) = i}0^a{x) . (4) 

Due to the periodicity in 9, general observablcs Og = 
{0[^-)g\) can be represented in a Fourier decomposition, 



Oi. 



(5) 
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This impHes that the dual observables Oi change un- 
der a center sensitive gauge transformation Uz into z'O/. 
Hence every moment Oi with I E Z and I mod Nc ^ 
has to vanish in the center symmetric phase as it is pro- 
portional to a sum over center elements z, 



c Ol mod N^,0 ■ 



(6) 



Note that Ce — e^'^'^(L) reflects the boundary condi- 
tions ([3]) and fits into the definition of dual observables. 
The observables Og can either be evaluated in QCD with 
anti-periodic quarks @, 0,11, Hi, or in QCD^ with quarks 
having 0-dependent boundary conditions. 

In summary, the moments Oi with I mod ^ are 
order parameters for the confinement phase transition in 
QCD if evaluated in both, QCD and QCDg. In particular 
the first moment 0\ is an order parameter for all iVc, 



(7) 



For example, the dual Polyakov loop in QCD is (L). 

QCD at imaginary chemical potential - In the present 
work we mainly concentrate on QCD^. Its generating 
functional is 



(8) 



with (j>g ~ {A, i/jg^Tpg, ...) and J = {Jji,fj, f], ...). The dots 
stand for theghost fields and composite hadronic fields, 
see e.g. [13: El- Here, S denotes the standard QCD ac- 
tion and includes a Dirac action with 0-depcndent quark 
fields. With (g]) wc have 



V'e {ilp + rn) ipg 



ij[ip + m-2TT^joO]^, (9) 



where ^ = + g^). The rhs of ([9]) is nothing but 
the Dirac action with an imaginary chemical potential 
H = 2Tri9/f3. If 9 takes one of the center values 9z, we 
can define ip^' = ipg with Uz as in ^ and anti-periodic 
■0. We conclude that center phases dz can be absorbed 
in center transformations of the gauge field, A A'^^ . 
The generating functional at vanishing current J has the 
Robcrgc- Weiss (RW) periodicity [13], see also [i. [Tsl. [l6j . 



Zg[0] = Z, 



9+1/N, 



[0] 



(10) 



Furthermore, observables Og in QCDg arc RW- 
symmetric for J — 0. Hence only the center-symmetric 
Fourier coefficients Onj arc non- vanishing. In turn, a 
non-vanishing current for the gauge field breaks the 
RW-symmetry, and leads to O; 7^ for I mod Nc 0. 
This leads us to a simple and easily accessible confine- 
ment order parameter in QCDg, the dual density: 



>J\ ■= 



= / dOf 



ng = i/3 



dOf 



[J], 
(11) 



with 4>j = {cj))j. The density ng is the derivative of the 
partition function w.r.t. the chemical potential 2Tr 9/ f3, 

ng[M = J d''x{^-/o^p)g = ^dg\nZg[J] . (12) 

On the rhs of pT|) we have integrated by parts and made 
use of Zo[J] = Zi[J]. The dual density is propor- 
tional to the first moment of the grand canonical poten- 
tial in the presence of a gauge field background ip. Hence, 
it grows like at high temperatures, as the integrated 
6'-dcpendence is expected to be leading order. 

The above analysis for the dual density extends to 
general observables Og[4>j]. They constitute observables 
in different theories distinguished by the boundary 
condition. Hence the order parameters O in ([7]) in the 
presence of a fixed background 4>j = only vanish if 
QCDg is in the center symmetric phase for all boundary 
conditions. However, the transition temperature in 

QCD, Teonf = Teonf (0 = 0), is & lowCr bouud for Teo„f(0), 

see Fig. [31 Thus, the dual phase transition temperatures 
Tconf are identical with the physical one, Tconf = ^conf- 

Analytic properties - Some interesting properties of 
dual order parameters can be accessed analytically. To 
that end it is convenient to study observables in terms 
of the quantum effective action r[(/)] = / J0 — lnZ0[J]. 
The current J is given by J = 5T[(f>\/ 5(j)^ and vanishes 
on the equations of motion with ^-dependent mean value 
4>e = 4>J=o- For example, the dual density pip turns into 



d9( 



(13) 



As discussed above, we have — 0. This follows 

from the RW-symmetry of T[(f)g] which is a consequence 
of that of Zg. For its direct proof it is sufficient to ex- 
amine r[(/)] for constant gauge field configurations Aq in 
the Cartan subalgebra. For example, for A^c = 3 we have 
Cartan fields (3gAo = (p3T^/2 + ipsT^/2 with Gell-Mann 
matrices t^,t^. The 6'-dcpendcncc of the effective action 
r[(/)] originates in sums of terms with fermionic Matsub- 
ara frequencies stemming from Dq + 2ttT9: 2TTT{n+l/2+ 
9 + f3gAo/ {2Tr)). Most of the ^-dependence can be reab- 
sorbed in a ^-dependent gauge field ^o(^)- For Nc = 3 
this leads to Cartan components ipa = ^fa — 3(27r)0 and 
ips ^ ips- \/S{2n)9 of ^0(6'). Now, the RW-symmetry is 
explicit in the Matsubara frequencies 



2TrT {n+l- + 

z 47r 



i = l,...,A,, (14) 



where the $i's are the eigenvalues of the matrix 
2PgAo{9): 9 9 + 9z is absorbed in a center gauge 
transformation of the ipiS as well as in a shift of the 
Matsubara sum. Under this combined transformation 
the (^i's are invariant and so is the effective action. 
In particular we conclude that any expansion scheme 
based on fixed field variables ipi is form-invariant under 
9 ~* 9 + 9z- Moreover, the observables Og[<f>g] are in- 
variant, and hence ©[^e] = 0. In turn, observables 0[(f>] 
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Figure 1: Functional flow for the effective action: The lines 
denote full field dependent propagators. Crosses denote the 
cut-off insertion dtR- 



with ^-independent gauge field background ip are order 
parameters for confinement as such a background expHc- 
itly breaks the RW-symmetry. In particular this includes 
0[(j)] with ip ^ (p = (pe^Q and Lp = 0. 

Simple observables 0[(j)] follow directly from the ver- 
tices r("^[0] in QCDg. This includes the dual density 
as well as the dual chiral condensate with Og[(j)j] = 
J d^x {'ipeipe) J for either (pj = 1^9=0 and (pj = 0. The 
first case with (p relates to the lattice computations in 
QCD of dual order parameters B, Bl ■ The latter choice 
has been used implicitly in [9|, ll7[. An even simpler 
observable is the dual quark mass parameter M with 
Meifp] ^ tr T';^^[(j)]{p = 0). The specific choice Meicpe] is 
directly related to the pion decay constant f-^ in QCDg. 
A further prominent example is the modified Polyakov 



loop variable Lg 



„2mB 



L, 



(15) 



with Cq = {Lq). Eq. ([TS]) is invariant under 9 ^ 6 + 6^ 
at fixed ip, and hence /C[0e] = 0. However, L[(p\ — L\0\ 
simply is the Polyakov loop variable introduced in p^.[l9j 
as an order parameter for confinement. 

The representation of the Polyakov loop in (jl5p leads to 
an interesting observation: in phase-quenched QCDg we 
are left with the explicit 6'-dependence in the Matsubara 
frequencies (fTi|) . Thus, any observable O in ([7]) obeys 



d0t 



'Og[0]L[p]=d[0]L[p], (16) 



for ^-independent gauge field background p and vanish- 
ing quark and mesonic backgrounds. In fully dynamical 
QCDg the factorisation only holds approximately. 

QCD with functional methods - Our numerical compu- 
tations are performed within the functional RG approach 
to the quantum effective action Ffc. Here, k is an infrared 
cut-off scale below which quantum fluctuations are sup- 
pressed. For /c — > we regain the full quantum effective 
action F. The effective action F^ obeys the Wetterich 
equation poj . 



1 



9tFfe[0] - -Tr 



1 



Rk 



-dtRk , 



(17) 



with t = Ink /A and cut-off functions that provide in- 
frared cut-offs for all fields, for reviews on gauge theories 
see [1^ [ll| . The diagrammatic representation of (|17p is 
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Figure 2: The pion decay constant fniT)/ fTr{0), the dual den- 
sity n(T)/n(oo), and the Polyakov loop L[ip]{T) as functions 
of temperature, Xl = drL, Xduai = drfi. 



provided in Fig. [T] The RG fiow for the dual density fol- 
lows from (|13p and is directly related to the flow of the 
effective action. This is an important property as the flow 
for Ffc is least sensitive to the approximations involved. It 
also guarantees the maximal disentanglement of the dif- 
ferent fleld sectors, see Fig.[TJ We conclude that the dual 
order parameters, in particular n, are dominated by the 
quark loop. Moreover, in the present approach with dy- 
namical mesonic degrees of freedom the quark propagator 
is in leading order only sensitive to the chiral properties, 
see e.g. [2l|, l22]- Consequently, the confinement temper- 
ature in QCD, derived from the dual density n, has to 
agree approximately with the chiral critical temperature. 
Thus the chiral and confinement phase transitions are 
necessarily closely related. This observation is sustained 
by our explicit computations, see below. 

Numerical results & summary - For our study of two 
fiavour QCD in the chiral limit, we solve the fiow equa- 
tion for the effective action F b y c ombining results for 
the Yan g-M ills part of QCD [H,!!^, as well as the mat- 
ter part [21]. [2^. [2^ . The two sectors are coupled by the 
dynamical quark-gluon interaction. This setting incor- 
porates the confining properties of QCD [31 via the full 
momentum dependence of gluon and ghost propagators 
[23! . The results for pure Yang-Mills agree quantitatively 
with the corresponding lattice results. In the matter sec- 
tor mesonic degrees of freedom are dynamically included 
im, m, US] • Such a treatment of the matter sector al- 
ready provides quantitatively reliable results for the me- 
son spectrum, see e.g. [13]. It has been also successfully 
implemented for the phase diagram of one fiavour QCD 
at finite chemical potential [2^ . 

In Fig. [5] the temperature dependence of two order pa- 
rameters for confinement are shown, namely the Polyakov 
loop variable L[0\ and the dual density n[(j)\. The 
crossover temperature Tconf is determined by the peaks in 
the respective T-derivatives xl and Xduai- Interestingly 
the factorisation (|16p works quantitatively for the dual 
density in the full theory: n[t/)]/n[0] and L[(p\ agree on 
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Figure 3: Chiral (T^) and confinement temperature (Tconf) as 
functions of temperature and boundary angle 9. The width 
oi XL = OtL is displayed as a shaded area. The dots indicate 
the endpoints of the confinement RW-phase transitions. 



the percent level. We have checked further order param- 
eters such as the dual pion decay constant. We find that 
the crossover temperatures extracted from the dual den- 
sity, the dual pion decay constant and the Polyakov loop 
agree within a few MeV: Tduai « Tconf ~ 178 MeV. This 
provides further support for the quantitative reliability 
of the present approximation. 

In Fig.[2]wc also show the pion decay constant It is 
proportional to the quark mass parameter Mg evaluated 
at and is an order parameter for the chiral phase 
transition. For T ^ 0, approaches 90 MeV. For T > 
181 MeV the pion decay constant tends to zero 
and chiral symmetry is restored. We observe a second 
order phase transition, and the critical exponents such 
as v signal the 0(4)-universality class. Most importantly, 
the chiral phase transition and the confinement crossover 
temperature agree at vanishing chemical potential. 

An evaluation of the dual chiral condensate and the 
dual quark mass parameter in QCDg for vanishing gauge 
field background ip = has been implicitly performed in 



and [17| respectively. Evaluated at both, and (p, we 
find the expected periodicity of M in ^ 0-1-1, and no 



RW-symmetry. For = 1/2 it can be shown analytically 
that it grows with T^^^ for large T. In turn, for 9 = 
and (p = (pg=o it agrees with and vanishes for large 
T, see Fig. [21 Below the chiral phase transition temper- 
ature the mass parameter A1 is a smooth function of 
9. However, a box- like behavior emerges above T^, see 
also [l3|- Details will be presented elsewhere. Here we 
simply note, that A^e[0] is an expansion parameter of 
the effective action which only signals chiral symmetry 
breaking for (f> = (f>g. 

In Fig. [3] we show the phase diagram of QCDg. The 
confinement and the chiral temperatures lie close to each 
other for all imaginary chemical potentials. Their value 
at 6* = 1/6 is the endpoint (Trw,6Irw) « (210MeV, 1/6) 
of the corresponding RW phase transitions shown as a 
vertical line at = 1/6 in Fig. [3l Our results c omp are 
well to the lattice results [ll| . In the PN JL-model [H] the 
lattice results have been reproduced by adjusting model 
parameters connected to an eight quark interaction p^ . 
In our approach to full QCDg coinciding temperatures re- 
sult from the interplay of quantum fluctuations and are 
not adjusted by hand. An estimate of the correspond- 
ing quantum fluctuation within a Polyakov-quark-meson 
model also leads to coinciding critical temperatures at 
real chemical potential These results suggest that 

the differences between Tconf and at both, real and 
imaginary chemical potential, are mainly due to mean 
field or large Nc approximations. The relevance of this 
observation for the quarkyonic phase proposed in [27| will 
be discussed elsewhere. 

In summary our study suggests that the confinement 
and chiral critical temperatures Tconf and are dynam- 
ically related and agree within the error bars. At present, 
we extend our work to real chemical potential. This may 
help to shed some light on the current debate concerning 
lattice simulations at finite chemical potential. 
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